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Abstract. We study the Josephson effect in alkali atomic gases within the two- 
mode approximation and show that there is a correspondence between the Bogoliubov 
description and the harmonic limit of the phase representation. We demonstrate that 
00 ' the quanta of the Josephson plasmon can be identified with the Bogoliubov excitations 

I of the two-site Bose fluid. We thus establish a mapping between the Bogoliubov 

approximation for the many-body theory and the linearized pendulum Hamiltonian. 
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PACS numbers: 03.75.Fi, 74.50. -fr, 05.30.Jp 
1. Introduction 



B 

o 
o 

The Josephson effect and vortices are at the heart of superfiuid phenomena studied 
^ . for decades in superconductors and hquid Hehum. The achievement of Bose-Einstein 
■ condensation in trapped alkah atomic gases fl]] opened the way to further extend 
our understanding of these macroscopic quantum phenomena in a new system: cold 
dilute bosonic gases. Predicted forty years ago for Cooper pair tunneling from 
one superconductor to another through an insulating junction, the Josephson effect is, 
mathematically, a result of restricting the number of states available in the one-particle 
Hilbert space to two orbitals into which N bosons are distributed. In real life this is 
usually done by weakly connecting two superfiuids. This generates a wealth of physical 
phenomena that are essentially due to collective oscillations between the two superfiuids. 
For trapped alkali gases the two superfiuids can be either two condensates of the same 
species spatially separated by a potential barrier (external Josephson effect) |^,|[ or can 
correspond to spatially overlapping condensates made of atoms in two different hyperfine 
states (internal Josephson effect) 
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The standard Josephson Hamiltonian is the two-mode version of the Bose- Hubbard 
model, 

H = -^{a^b + b^a) + ^ [{a^af + (b^b)'] , (1) 

with the Josephson couphng energy Ej and the charging energy Ec- Intrinsic to the two- 
mode approximation is the assumption that the parameters Ej and Ec are constant, i.e. 
independent of the number of particles in a given state. Thus we assume implicitly that 
fluctuations in the relative particle number do not change the shape of the condensate 
wavefunctions significantly. 

In the case of identical wells, a and b are the annihilation operators corresponding to 
the condensate wavefunctions of the left and right well; an evaluation of the parameters 
Ej and Ec has been given in Ref. 0. In the case of the internal Josephson effect, atoms 
in two different hyperfine states |1) and |2) are confined in an optical or magnetic trap; 
to implement the Josephson coupling term a Raman transition is driven between the 
two levels ^ . The many-body Hamiltonian describing the mixture of the two quantum 
fluids with scattering lengths an = can be reduced to the effective form (|lD with 
parameters Ec oc (an — au) and Ej = Nhfl/2 0, where Q is the Rabi frequency 
and the detuning is assumed to be zero. 

Coherent particle exchange between two different hyperfine states has been observed 
101 in the Rabi limit [Q of the internal Josephson effect. Then the two modes a and 
b can be for example the hyperfine states \E=2,mF=l) and |-F=l,m^=— 1) of ^''Rb 
for which the robustness of phase coherence has been already experimentally checked 
out in magnetic traps p[ using a combination of microwave and rf fields to drive the 
two-photon transition, or F = 1, mp = ±1 states of ^^Na which are miscible and can 
be trapped in optical traps 0. 

A vertical array of cold atoms trapped in the anti-nodes of an optical standing wave 
was used to create an analog of the ac Josephson effect and a lattice of quasi-one- 
dimensional confining tubes formed by a pair of laser fields was recently employed to 
investigate the phase coherence between neighboring sites 0. The physics of these 
experiments is based on the macroscopic coherent tunneling of atoms between the 
wells, and is very similar to that of the Josephson effect between two superconductors 
connected by an insulating junction |jTl[ . 



The model Hamiltonian (|ID has been studied in the classical and the quantum 
regime by using different mathematical techniques. Most popular are the phase 
representation |jl2[, the angular momentum representation [|13[, and the time-dependent 
Gross- P it aevskii (GP) equation in either its reduced phase-number formulation [p|,p!4| 
or in its hydrodynamic version |15|. As expected on very general grounds [p!6| ,p!7[] and 
also checked numerically |T8|, both GP approaches yield the same excitation spectrum 
as the one obtained by solving directly the time-independent Bogoliubov-deGennes 
equations. Up to now, only the relationship between the phase representation and 



the angular momentum representations has been pointed out [p!3|, O, |T9[ . In this 



paper we focus instead on the equivalence between the harmonic limit of the phase 
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representation and the Bogoliubov approximation |^0[ of the Hamiltonian (|1]). In 
particular, we estabhsh the relationship between the creation and annihilation operators 
in the phase representation and the Bogoliubov operators, and show that the quanta of 
small oscillations of the Josephson pendulum are Bogoliubov excitations in disguise. 

2. Phase representation and Josephson plasmon 

Let us introduce the relative number operator 



1 



n = -[a'a 



2."^- b'b), (2) 

and let N = a'' a + 6^6 be the total number of atoms which will be treated as a c- number 
since it is conserved. Then a simple way to obtain the phase representation of the 
Josephson Hamiltonian ([l|) is to insert a polar decomposition of the operators a and b, 

a = ^N/2 + ne-''P'\ h = ^N/2-ne''^'\ (3) 

This yields the momentum-shortened pendulum Hamiltonian 0,|T1| whose classical limit 
reads 



1 / 4^2 1 

H = -E,N^ - Ej^l -—COSV + -E,n\ (4) 

A quantum version can be obtained from symmetrizing the second term. This procedure, 
initiated by Dirac for the one-mode field operators in his original description of the 
electromagnetic field, is not without problems when it comes to finding explicitly a 
Hermitian operator tp [^]. A similar idea is often invoked in condensed matter physics: 



It is argued that, in certain conditions and generally for systems with large numbers of 
particles, the U{1) gauge symmetry is spontaneously broken so that the mode operators 
a and h can be replaced with well-defined complex numbers in an amplitude-phase 
representation, yielding the form (^). 

In the following we will adopt a definition of the relative phase operator inspired 



by Ref. p2[. We consider the states 

1 



1 N/2 

^v) = ^WT^ E ^^|n), (5) 



-Nj2 

where \n) is a shorthand for \N/2 + n^N/2 — n) and the phase = 2ttp/{N + 1), 
p = —N/2, . . . , N/2 has a discrete structure. Then we can define a phase operator by 

N/2 N/2-1 

e^^= Y: e^^''\pp){pp\= Yl \n + l){n\ + \-N/2){N/2\. (6) 

p=-N/2 n=-N/2 

It is easy to check that the Lerner criterion [n, e**^] = e*'^ is satisfied so one can also write 
[(p,n] = i. The eigenstates \ipp) satisfy the orthogonahty relation {ipp\ipp') = 5ppi and 
form a complete set, Yl,p=-N/2 \Vp){^p\ — ^i which makes them suitable for the definition 
of a phase representation. This construction solves the problem of defining a (relative) 
phase operator by imposing a discrete spectrum for it. In the limit of large the phase 
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spectrum becomes quasi-continuus and one may define a derivative d/d(f to obtain for 
tfie number operator tfie pfiase representation n = —id /dip. 

The representation given by tlie pfiase states \ipp) allows to write the Hamiltonian 
(|^) in the form of the momentum-shortened pendulum Hamiltonian (^). In the limit of 
small oscillations, this effective Hamiltonian becomes 

H ^ ^E,N^ -Ej + ^Ej^' + ^E,n^ (7) 

with the effective charging energy 

AEi 

Ec^E,+ ^-§^. (8) 
Equation (|^) describes a harmonic oscillator with the frequency 



000 = yEjEjh (9) 
and the root mean square of the number and phase fluctuations of its ground state 

fulfill the minimum uncertainty relation AnAip = 1/2. Writing n and (p in terms of 
creation and annihilation operators, 

n = (a + at) j^-J. j , (11) 

^ = ^{a-a^)[^^J , (12) 

where a and at fulfill the bosonic commutation relation [a, at] = 1, diagonalizes the 
Hamiltonian, i.e. brings it to the form 

H = ^E,N^ -Ej + huo (ata + ^) . (13) 

It is useful to define three regimes [^[1^ for the Josephson two- mode Hamiltonian 
according to the interaction strength E^, namely the Rabi regime E^ <^ Ej/N"^, the 
Josephson regime E^/N"^ <^ E^. <^ Ej, and the Fock regime Ej <C Ec. In the Rabi 
regime, the atoms are all in the bonding state, but behave independently. The phase is 
well-defined and the excitation is the promotion of a single atom to the anti-bondig state. 
In the Josephson regime, the grond state has still a well-defined phase, but the excitation 
forms a collective motion, the Josephson plasmon with the plasma frequency \/EjE^/h. 
In the Fock regime, the Josephson link is dominated by the interaction energy and n is 
a good quantum number. Therefore the ground state has a well-defined atom number 
on each side, the phase is completely undefined, and the harmonic approximation (^ is 
no longer appropriate. 

We close this section with a comparision between this and similar phase represen- 
tations that can be found in the literature. In |^ the phase is allowed to take continuus 



values from the beginning. In |^ the problems of defining a Hermitian phase operator 
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are overcome by introducing an over-complete phase representation, in a basis of phase- 
coherent states. Far from the Rabi regime, this representation becomes equivalent with 
ours (cf. Appendix A). Finally, in [2^] a number representation is used, which obeys 
the same commutator relations but with "position" and "momentum" interchanged (cf. 
Appendix B). 

3. The Bogoliubov approximation 



To implement the Bogolioubov approximation |^ for the Hamiltonian (|l|), we first 
determine the condensate wave function. As an ansatz we use the most general many- 
body state for atoms all occupying the same mode. 



1^) 



cost 



N 



I vac). 



Its energy energy expectation value is 



{'^\H\m) = -Ej sm2e cos (p + -E^N^ (2 - sin^ 2^ 



(14) 



(15) 



which becomes minimal for 6 = 7i/4 and 
bonding state, so that the operator 
a + b 



Co 



V2 



[co, cl] 



0. Thus, the condensation occurs in the 

(16) 



destroys a condensate atom. The remaing mode 
a — b 



Cl 



V2 



[17) 



is orthogonal to cq and consequently the commutation relations [q, Cj] = Sij hold true. 

The central assumption of the Bogoliubov approximation is that one can replace 
the operator cq by (A^ — cjci)^/^ - 



A^ 



CiCi/yN keeping only terms up to second 
order in ci. This results in the Hamiltonian 

H = ^E,N^ -Ej+ Qe.AT + 2^) c\c, + ^E,N [c,c^ + c\c\) . (18) 

We now employ the symplectic transformation 



Cl 
Cl 



wy — f 7 



■u 7' — 1> 7, 

where the ansatz u = coshx, v = sinhx ensures |up 
commutation relation [7,7^] = 1. The choice 

E. 



(19) 
(20) 

1 and, thus, the canonical 



tanh 2x 



E, + 8Ej/m 



brings the Hamiltonian ([T8|) to the form 

H = ^E,N{N -l)-Ej(l + ^^+ hcoo (7^7 + I) ■ 



(21) 



(22) 
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The ground state energy of this Hamiltonian is shgthly lower than the mean-field energy 
EcN"^ /8—Ej; this refiects the role of interactions, which in general distribute particles on 
modes other than the condensation state. In the limit ^ 1 this ground state energy 
becomes the same as that predicted by the phase representation. Also, the diagonal 
form of the Hamiltonian indicates that the energy of the Bogoliubov quasiparticles is 
the same as that of the Josephson- Rabi oscillator in eq. (|^) . 

The corresponding Bogoliubov ground state is defined by 7|BdG) = 0. Its depletion 
number, i. e. the number of atoms that does not reside in the one-particle ground state, 
A''{ = (cjci) is easily evaluated to read 

11/^. 1 y 




NAif - —= . (23) 



A^i' = sinh^ Y = — A 1 \ I — = - I V i V l^ijj = i 

' ^ 8\Ej 2N^ 2 A\ ^ VNAif) 

The condition for the applicability of the Bogoliubov approximation is A'^^ <^ A^, which 

means that it is not valid in the regime Ej < Ec- Comparing with Eq. (|10D yields that 

the Bogoliubov approximation breaks down when the phase is not well defined. In this 

case the depletion is so large that the Penrose-Osanger criterion is not satisfied anymore; 

the notion of a single condensate is no longer appropriate and the ground state of the 

system will be fragmented. Indeed, we find that the 2x2 one-particle density matrix has 

elements {a'^a) = (b'^b) = N/2 and (a^b) = (b'^a) = N/2 — N[ so to have condensation 

on the state {a + b)/ ^J2 one needs to make sure that the off-diagonal elements in the 

one-particle density matrix are of the order iV, i.e. the depletion number iV( is negligible 

with respect to A^. 

The structure of the Bogoliubov ground state |BdG) can be obtained by noticing 
that the transformation (|l9l),(p0|) is a squeezing transformation ^ 



7 = ci cosh X + 4 sinh x = S{x)c\S\x) (24) 
with the squeezing operator 

Six) = exp [^x{cl - )] . (25) 

The ground state structure is then 

|BdG)=5(x)|GP), (26) 

where |GP) denotes the Gross-Pitaevskii ground state. 

Let us establish a relation between the corresponding creation and annihilation 

operators by writing the number difference operator n in terms of the Bogoliubov 
operators, 

n = l(ata - b^b) = ^(cjco + cjci) ^ i= (^|^) ^ (7 + 7^). (27) 

To obtain the final expression, we have again used (cq) ~ and e~^^ = AEj /N'^E^. 
The operator conjugate to n is uniquely defined by the commutation relation [yj, rt\ = i 
and must therefore read 

# (7-V)- (28) 
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This demonstrates that the operator sets and 7,7^ are identicaL Therefore 
the Hamiltonian ( ]T3| ) of the hnearized pendulum is identical with the Bogoliubov 
Hamiltonian (^) and, consequently, the Josephson plasmon can be viewed at as a 
Bogoliubov quasiparticle. 



4. Concluding remarks 

In the Rabi and in the Josephson regime, a split condensate has a well-defined relative 
phase and that is why Bogoliubov theory works. With increasing effective interaction Ec, 
the uncertainty in the relative particle number decreases from /2 in the Rabi regime 
via [Ej / AEcY^^ in the Josephson regime to a value much smaller than unity in the Fock 
regime. At the same time, the fluctuations of the relative phase keep growing until the 
phase becomes completely undefined. In the Bogoliubov approach, this corresponds to 
an increasingly larger depletion that finally becomes of order and, thus, violates the 
condition that most of the atoms have to reside in the same one-particle state. 

In conclusion, we have established the equivalence between the harmonic limit of 
the phase representation and the Bogoliubov approximation in both, the Rabi and the 
Josephson regimes. The quanta of the Josephson-Rabi oscillator are the quasiparticles 
of the Bogoliubov theory. The vacuum of the Bogoliubov theory is the ground state 
of the phase-number harmonic oscillator. Thus, despite their different mathematical 
appearance, both approaches can be mapped onto each other and describe the same 
physics. 



Appendix A. The phase-coherent states representation 

The states (|^) are not the only possible meaningful phase states. Another widely used 
choice is the over-complete set 

1^^) = -jL= \a^e-'"^ + fote^^/^l"^ |vac), (A.l) 

where 6' is a continuus variable. In this representation and for Ej <^ N'^E^. the 
Hamiltonian reads 

H = ^E,N' - E J cose - IEc^, (A.2) 

which agrees with the pendulum Hamiltonian in the Josephson regime and n -C A^. 
The reason for this agreement is that for n <^ N the coefficients of the expansion of 
( |A.1| ) in the basis \n) can be approximated by 

|(r2|^)| oce-"'/^. (A.3) 

Since in the Josephson regime the number fluctuations are much smaller than V^, the 
Hilbert space is explored for n only in the range between ±^/N, and the coefficients 
( |A.3| ) become fiat. But this is precisely what characterizes the phase states (|]), so this 
argument proves that indeed in this regime the two descriptions become identical. 
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Appendix B. The number representation 

Instead of working in a basis of pliase states one can as well decompose into tlie 
number states \n), 

N/2 
n=-N/2 

and assume that \I/(n) changes smoothly between consecutive values of n. Going from 
one representation to another is achieved by the Fourier sum 

for wavefunctions, while for phase operator acts as a derivative, (f = id/dn. This 
brings the Hamiltonian to the form but now with ip being a derivative. After 
a linearization, we obtain a harmonic oscillator where "position" and "momentum" are 
interchanged. 

Another route, followed in [24|, is to start directly from the two-mode Hamiltonian 
(|I|) and to decompose it into the number states |n) to obtain 



With the operator identity 



cos li-^] ^(n) = - [^(n + 1) + ^(n - 1)] (B.4) 
\ an I 2 



follows the desired expression. 
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